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1 Introduction 

Different cosmological and astrophysical observations point towards the existence of Dark 
Matter (DM), a non-electromagnetic interacting form of matter that is responsible for struc¬ 
ture formation and other gravity-related phenomena. Among these, we may cite the observed 
high velocity dispersion in clusters of galaxies, the confinement of the hot gas detected in 
these objects and the flat rotation curve of spiral galaxies. In particular, the modelling of 
the rotation curve of the Milky Way requires, besides the known baryonic components, the 
presence of a massive dark halo [1, 2]. 

The nature of DM is one of the greatest challenges in cosmology and physics today. Since 
standard model relics are not able to explain the required cosmic abundance, DM is possibly 
related to extensions or modifications of the canonical model or can even require new physics. 
The most promising particle candidates are WIMPs (Weakly Interacting Massive Particles), 
which are interesting because a stable relic with mass close to the electroweak scale (~ 100 
GeV) has an expected abundance today able to explain the observations (i.e. Ddm ~ 0.3). 
Unfortunately, despite the great efforts deployed by different experimental groups trying to 
detect signals from the interaction between these particles and those of the standard model, 
until now no positive result has been reported. For a review on DM particles, see for instance 
ref. [3]. 

In order to trigger structure formation, DM must behave like a pressureless fluid or, 
in other words, to have a quite small velocity dispersion, reason why it is called cold DM 
(CDM). However, in the early universe WIMPs are still coupled to the primordial plasma, 
having thus a non-negligible pressure. This is consequence of interactions with standard 
model particles (essentially leptons) that produce damping effects in the primordial density 
fluctuation power spectrum on scales smaller than a damping (or diffusion) scale Ad, which 
depends on the interaction processes which characterise the DM particle model. 

For example, a DM-nentrino conpling would erase small DM structures due to colli- 
sional damping effects and late neutrino deconpling, see refs. [4-6]. Simulations of non-linear 
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structure formation in such models were analysed in ref. [7]. In refs. [8, 9], the authors showed 
that a DM-photon or DM-neutrino coupling, with a cross section compatible with the CMB, 
could alleviate the missing satellite problem [10]. Moreover, effects related to velocity dis¬ 
persion have proved to be relevant, e.g. in the Tseliakhovich-Hirata effect [11], where the 
relative motion among baryons and DM is considered. Se also refs. [12-15]. 

After kinetic decoupling, WIMPs “free stream” but still possess a non-negligible velocity 
dispersion that fixes the free-streaming scale Afs below which the gravitational collapse is 
suppressed. In reference [16] the authors derive analytically an approximate expression for 
the free-streaming scale, namely, Afs ~ 1 pc and for the diffusion scale Ad ~ 10“^ pc. A more 
detailed analysis can be found in [17]. A different approach is adopted in ref. [18] where the 
authors consider first order fluctuations in the velocity dispersion, after kinetic decoupling. 
These may behave as pressure perturbations, allowing one to define a Jeans scale for WIMPs 
that is slightly larger than Afg. Consequently, the cut-off of the matter fluctuation power 
spectrum at small scales is controlled by the Jeans scale deriving from velocity dispersion 
fluctuations. 

If DM particles decouple relativistically from the cosmic plasma, they constitute the 
model dubbed Warm Dark Matter (WDM) [19-26]. A thermal origin implies that particle 
masses are of the order of few eV in order to explain the observed abundance while a non- 
thermal origin allows masses of the order of few keV [27-30]. See also refs. [31-33]. 

WDM may alleviate some problems that plague CDM at small scales such as the missing 
satellite problem [10], the core-cusp problem [34] or the too big to fail problem [35]. WDM 
effects related to the the core-cusp problem were studied by high resolution simulations per¬ 
formed in ref. [36]. These simulations have shown that the observed size of cores can be 
reproduced if the mass of WDM particles is about 0.1 keV. However such a low mass pre¬ 
vents the formation of dwarf galaxies since introduces a cut-off in the matter power spectrum 
at scales less than 10^® A7 q. 

On one hand, supporters of WDM, see for example refs. [22-24], argue that quantum 
effects on scales smaller that 100 pc are relevant, and when properly taken into account, e.g. 
via a Thomas-Fermi approach, permit to predict sizes of galactic cores, masses of galaxies, 
velocity dispersions and density profiles in agreement with observations. On the other hand, 
recent investigations of WDM based on the flatness problem of the HI velocity-width function 
in the local universe and Lyman-a forest [25], fail to alleviate small scale problems related 
to structure formation. According to [25], this failure is probably due to the cut-off in the 
linear power spectrum that is too steep to simultaneously reproduce the Lyman-a data and 
to match dwarf galaxy properties. See also the very recent work [37] for the tightest bound 
to date on WDM mass. 

In this paper, the effects of velocity dispersion on the linear growth of DM fluctuations 
are analysed. We start from the Vlasov-Einstein equation for DM, building a hierarchy of 
equations that are truncated at the second order momentum, by neglecting contributions 
of the order of 0{p^/mfl), where p and m are respectively the proper momentum and the 
mass of a DM particle. Then, a modified GAME code [38] was used to compute the linear 
matter power spectrum. The derived results are compared with the matter power spectrum 
derived from a WDM model constituted by particles of mass ~ 3.3 keV. Both models are 
able to produce a cut-off in scales of about ~ lOh Mpc“^ but for different physical reasons. 
We speculate if velocity dispersion effects would be more appropriate for solving small-scale 
issues present in the DM paradigm. In this case, a late non-thermal DM origin is required 
and our picture becomes similar to the Late-Forming scenario, in which DM particles are 
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produced by the decay of the kinetic energy of a scalar held in a metastable vacnum. 

The paper is organized as follows. In sec. 2 the equations describing the evolution of 
small CDM huctuations in the presence of a non-zero, small velocity dispersion are presented. 
In sec. 3 the modihcations of the DM power spectrum and the CMB TT power spectrum 
induced by the presence of velocity dispersion huctuations are discussed. Finally, in sec 4 
the main results and conclusions are given. 


2 Evolution equations for small dark matter fluctuations 


In this section we derive the equations governing the evolution of DM huctuations in presence 
of a velocity dispersion. We write the perturbed Friedmann-Lemaitre-Robertson-Walker 
metric in the following way: 

ds^ = a^(r) [—-|- {5ij -|- hij) dx^'dx^] , (2.1) 


where h^j, focusing only on scalar perturbations, is written in the following form as an inverse 
Fourier transform: 


hjj(x, r) = J d^fcexp(ik-x) 


kikjh{k,T) + ( kikj - -6ij 


6rj(k,T) 


( 2 . 2 ) 


where x is the position 3-vector, k is the wavenumber 3-vector with modulus k and direction 
ki, and h and r] are the scalar perturbations of the metric, i.e. the trace and traceless part, 
respectively. Note that, by dehnition, 5ijk^k^ = 1. 

We are going to determine the evolution equation for DM with velocity dispersion by 
taking up to the second momentum of the Boltzmann equation. Dehning / = /o + /i as 
the one-particle distribution function, separated in its background unperturbed value and a 
small huctuation about it, the perturbed Boltzmann equation can be computed, see ref. [39], 
as follows: 


dr 


a d/i p M r _L ^-^0 

-Pir~ + F • n /i + 
a op E op 


^(k-n)^ 


= 0 , 


(2.3) 


where p® = pi is the proper momentum; h® is the unit vector describing the direction of p®, 
i.e. p® = ph®, where p is the modulus of p®, i.e. p^ = djjp®p^ (remember that the definition 
of proper momentum incorporates the metric); E'^ = p‘^ + m? is the energy; k is the wave- 
vector coming from Fourier transformation and the dot denotes derivation with respect to 
the conformal time r. 

In the following, we consider momenta of Boltzmann equation (2.3) in order to derive 
the evolution equations for the perturbative quantities. This task is similar to the usual 
one for photons and neutrino. On the other hand, in our model for DM we will consider 
non-relativistic “quasi-cold” particles, in the sense that they have a non-vanishing bnt small 
velocity dispersion. 

We implement this condition by employing the following expansion for the particle 
energy: 


E = m 



2m? 


+ 0{p^/m!^) , 


(2.4) 


and by neglecting terms of order O (p^/m^). We will use the symbol ~ when making this 
truncation. This will allow us to truncate the hierarchy of the Vlasov equation at the third 
order momentum. 
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2.1 Zero-order momentum of the perturbed Boltzmauu equatiou 

Taking the zero-order momentum of eq. (2.3) means to multiply it by d?p and to integrate 
in the proper momentum space. The result is: 


drii a 1 i h + 6r] 

— -h 3-ni -I- ikiTiov - 3nor] -\ -- —uq = 0 , 

OT a 2 

where we used the definitions 

n = no + ni , ^ j / ^^P^^ ’ 

for the particle number density^ and 


nov 




(2.5) 


( 2 . 6 ) 


(2.7) 


for the bulk velocity, which is a pure first-order quantity. We also used the geometrical 
property 

J dnn^h^ = , ( 2 . 8 ) 

in order to integrate by parts the last term in eq. (2.3). Defining the number density contrast 
as (5n = ni/no, we can rewrite eq. (2.5) as: 


6n T ikiV^ T 77 — 0 , 


(2.9) 


where we also used the background result ho-|-3dno/a = 0, i.e. particle number conservation. 
For CDM, 5n = S = pi/po, but in our case, considering a non-vanishing velocity dispersion, 
that is no longer true. The energy density can be calculated as follows: 

p = J dP’pEf = I d^p^yp^ + m^{fo + fi) J d^p (/o -k /i) . (2.10) 

Using eq. (2.6), the above integrals can be rewritten as follows: 

p = PQ + pi, po = nom + ^noma'^ , pi = nim + ^nomvf , (2.11) 

where cjg is the background velocity dispersion and vf its perturbative counterpart, defined 
as follows: 

nocjg = j d^p ^2/o ~ j d^P^fo , riQvl = J d^P^fi ~ j d^P^fi ■ (2.12) 

From eq. (2.11) we can read out the density contrast: 


^ _ Pi _ nim + mnovjj2 _5n + v\j2 
Pq nom + mnQal/2 l-|-cr^/2 

and rewrite eq. (2.9) as follows: 


1 + 




d - -al5 - + ikiE + - = 0 , 


2 dr 


h 

2 


(2.13) 


(2.14) 


where we used the result cTq oc a 


'^Note that the number of spin states, usually denoted by g, is incorporated in /. 
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2.2 First-order momentum of the perturbed Boltzmauu equatiou 

In order to compute the first-order momentum of eq. (2.3), we multiply it by d^ppn''/E and 
integrate. This procedure gives the following result, using also eq. (2.7): 

^ + ikj j d^p^fdn^ fi = 0 , (2.15) 


where we have used the fact that the integration of /o with an odd number of unit vectors h* 
is vanishing because of the isotropy of the background solution. The last integral of eq. (2.15) 
defines the perturbed velocity dispersion tensor: 

/ 2 r 2 

d^p^fdn^fi ~ J d^p^fdWfi , (2.16) 

consistent with eq. (2.12), since dijfdn^ = 1. Contracting eq. (2.15) with ki and using the 
background solution for no, one obtains: 


kiE + -fcjU* -|- ikikj{v‘fy^ = 0 . (2-17) 

2.3 Secoud-order momentum of the perturbed Boltzmann equation 

In order to compute the second-order momentum of eq. (2.3), we multiply it by p^rdWd^p/ 
and integrate. One hnds: 


(9[no(uf 




dr 


+ f)-no{viy^ - 5r]no{cro) 




d^p^foEnPh^fi^ = 0 , (2.18) 


where we used the dehnition 

noic^oT^ = J d^p^fofdh^ ~ J d^p^foEE , (2.19) 

for the background velocity dispersion tensor. In order to calculate the last integral of 
eq. (2.18), we need to compute first the following integral: 


pjim ^ j _ ( 2 . 20 ) 

From eq. (2.8) it is not difficult to determine that 

jijim ^ ^ . ( 2 . 21 ) 

So we can now cast eq. (2.18) in the following form: 

- 5y{a^y^ + + 2EkA a^ = 0. (2.22) 

dr a b V / 

The Boltzmann hierarchy is thus truncated and we have a finite number of equation (4 instead 
of the usual 2 for CDM) describing our “quasi-cold” DM model. 
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2.4 Summary of the equations 

Considering scalar perturbations only for the velocity, i.e. ikiV^ = kv, we can summarize the 
three evolution equations (2.9), (2.17) and (2.22) found so far: 


6n kv — 0 , 


+ -V — kkikjivly^ = 0 


a 


d{vi 


2^3 \ 


dr a D 


5*^' + 2FF) (jg = 0 


Now, we expand (nf)*-^ in its trace and traceless part, i.e. 

[vlf = ^ 

so that eq. (2.24) can be cast in the following form: 

« ^2 T i-r 

V -\— V — vf — fcn = 0 , 
CL 3 

where we have defined the anisotropic shear as 


n = kikjT.^^ . 


(2.23) 

(2.24) 

(2.25) 

(2.26) 

(2.27) 

(2.28) 


Moreover, we can separate eq. (2.25) in its trace and traceless contributions. The former is: 

d{vl) 


Qj 2 2 r\ 

+ 2-'^! + "^'^0 = 0 ! 
dr a 6 


whereas the traceless part is: 


h + 6fj 


+ 2-S*^ + 
a 3 


-CTn 


) = 0 
o 


(2.29) 


(2.30) 


Contracting the above equation with kikj we obtain an evolution equation for H: 


■ 2 h + 6f] r, 


2 


a 

a 


6 


h 

+ 2 

d k 2 , 

- -V - -vi - k 

a 3 

2^ - 5h 


a 


-'1 


6 


d(vi 

dr 


n + 2“n + lhtl’'-^ 

a 3 3 


(2.31) 


vi,Il,h and r]: 


= 0, 

(2.32) 

= 0, 

(2.33) 

= 0, 

(2.34) 

i = 0. 

(2.35) 
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This system has to be coupled with the equations for the other components and with Einstein 
equations. Based on our assumptions of collision-less Boltzmann equation and truncation of 
the particle energy terms of order 0{p^/m^), they describe the evolution of non-relativistic 
DM after kinetic decoupling. Erom Ref. [39], the Einstein equations can be cast in the 
following form: 

--% = AnGa^dT'^o , 

2 o 

(2.36) 

k^fl = A7rGa^ik^5T^i , 

(2.37) 

h + 2 — h — 2k^ri = —87rGa^(5T% , 
a 

(2.38) 

h + &fi + 2^{h + 677 ) - 2k^r] = 2A7rGa‘^kiy{5T^j - 5T\6^j/3) , 

(2.39) 

where the perturbed stress energy tensor introduced above is the total one. 
part is changed in our treatment, in the following way: 

Only the DM 

<5r%(DM) = - J d?pEh = -Pi , 

(2.40) 

(5r°j(DM) = J d^pphifi ~ nomvi , 

(2.41) 

5Tj(DM) - j d^P ^ ~ nomivjYj . 

(2.42) 


Therefore, the perturbed DM stress energy tensor combinations which enter the Einstein 
equations on the right hand sides are: 


6T%{BM) = -po6 , (2.43) 

ik^6T\{BM) = Po—^ , (2.44) 

6T\{DM) = , (2.45) 

i + (Tq / z 

ky{STj - 5r^5V3)(DM) = po—^^ . (2.46) 

1 + a ^/2 


Combining eqs. (2.23), (2.33) and (2.34) it is possible to find the following second-order 
equation for 6n- 


d • 2 

H-T 

a 


TT+ n 


h ah 

+ 7T 4—7T ~ 0 • 
2 a 2 


(2.47) 


The last two terms on the Ihs are the usual gravitational potential contribution, which from 
eqs. (2.36) and (2.38) reads: 


1 + - | = AttGo^ {6T% - 6T\) . (2.48) 

2 a 2 

Note the term multiplying in eq. (2.47). In general, it can be cast as a sum of a con¬ 
tribution dependent from Sn, thereby acting as an adiabatic pressure perturbation, plus an 
entropy contribution. Both affect the evolution of 6n on small scales in a manner we show 
in the next section. The evolution equation for 5 can be found from eq. (2.47) using the 
relation (2.13). 
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2.5 Fluid interpretation 

In Ref. [39] the following two equations, coming from the stress-energy conservation equation 
can be found: 


<5 = 

e = 


-(1 + w) 




—-(1 — 2 > w )0 
a 


w 

1 + w 


e + 


-r-tc 1 0 

5p 

6P/6p 
1 -I- rc ' 


k^6 - k^a . 


(2.49) 

(2.50) 


These equations are equivalent to the first two momenta of the Boltzmann equation for 
DM that we have considered previously, viz. eqs. (2.32) and (2.33), provided the following 
identifications are made: kv = {1 + 5f7Q/6)0 and IT = —(1 + 5f7Q/6)cr. The equation of state 
parameter and the pressure perturbation are given as follows as functions of cJq and 


Pq _ q-g/3 5P _ uf/3 

po 1 + ( 7^/2 ’ po I + all2 ■ 


(2.51) 


Notice how velocity dispersion endows the fluid with a time-dependent equation of state and 
how vl plays the role of a pressure perturbation. The effective speed of sound = Pi/pi 
and the adiabatic speed of sound = Pd/ po have the following forms: 


CefT = 


vl/2, 


5aoV3 

5n + vl/2 ’ “ 3 + hal/2 

When perturbations are adiabatic, then = c^, which gives: 

- 3^0<^n - 3^0^^ 5^2/6> 


(2.52) 


(2.53) 


where we used eq. (2.13). With this condition, eq. (2.34) becomes: 


<5n + — — 0 , (2-54) 

which, when compared to eq. (2.9), implies iPvi = kv = 0. In turn, this implies from 
eqs. (2.33) and (2.35) that 

y + n = 0 , (2.55) 

al{3h + 8p) = 0. (2.56) 


The last equation implies either a vanishing velocity dispersion, thereby frustrating our pur¬ 
pose, or it gives an unnatural relation between the geometric perturbations, i.e. 3h = —8?). 
Therefore, the conclusion is that the perturbations cannot be adiabatic. We investigate this 
issue in some detail in the next subsection. 


2.6 Adiabaticity 

First of all, we write down the first law of thermodynamics: 


TdS = dU + PdV , 


(2.57) 
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where the variables have the usual meaning. Now consider U = pV and V = N/n, where 
N is the particle number, which we consider hxed i.e. no DM particle are being created (or 
annihilated). With these positions we can write: 


Tds = ^- , (2.58) 

n 

where s = S/N is the entropy per particle. Remember that we assumed conserved particle 
number, therefore N could enter in the differential dS. The density, pressure and number 
density which appear in the above equation are the total ones. If we ask adiabaticity for each 
single component, one gets for DM: 


5 

^ + P/P 


(2.59) 


and using eq. (2.51) one obtains again eqs. (2.13) and (2.53), along with the related problems 
we already mentioned. 

If we demand that the total entropy fluctuation is vanishing, we get that: 


TT^ “ 4*^^ ’ 


(2.60) 


and similar equations for the other components (neutrinos and baryons). Now, on the back¬ 
ground, energy density and particle number conservation ensures that TqSq = 0, i.e. adia¬ 
baticity. Note that TqIiti = a^j^. At the hrst perturbative order, using also the equations 
for the perturbations that we derived above, one can find that 

—Si = . (2.61) 

m 6 

Therefore, the initial condition n = 0 also implies that the initial time-derivative of DM 
entropy is vanishing. 


2.7 Initial conditions 


In order to properly modify CAMB we need not only to include the new equations contain¬ 
ing the velocity dispersion, but also to treat properly the initial conditions (IC). From the 
discussion of the previous subsections, in our model we cannot have adiabatic IC. 

Following e.g. [39] and [40] it is not difficult to see that the IC for h and p are not 
modified by the presence of a velocity dispersion since we are deep in the era where photons 
and neutrinos dominate and in our model we leave these components untouched. Therefore, 
d-y = = —‘Ihj'i. 

Using eqs. (2.32)-(2.35) in the limit /cr —>■ 0, we obtain the following IC for DM with 
velocity dispersion: 


1-^ 


A 2 ^ 

'*=-” 1-2 


1 + 


5a, 


6 

n = 0 , 

2 5/12 

"^ = -18^°’ 

TT _ 2h + 6p^2 

li —- (7r\ . 

9 3 ° 


Of course, for a vanishing Uq one recovers the usual adiabatic IC. 


(2.62) 

(2.63) 

(2.64) 

(2.65) 
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3 Numerical results 


We use as free parameter of the modified equations in CAMB the velocity dispersion of DM 
at the present moment, which we dub 0^. Since the backgronnd velocity dispersion scales as 



^0(kd)«M 


(3.1) 


where and Okd are respectively the velocity dispersion and the scale factor at the kinetic 

decoupling, is dehned as 

“ '^0(kd)®kd • (3-2) 

If DM has been thermally produced, then our approximation of non-relativistic particles 
allows us to use a Maxwell-Boltzmann distribution for DM particles and to write = 

3Tkd/m, where T^d is the cosmic plasma temperature at kinetic decoupling and m is the 
particle mass. In this case we have for /3^: 


,2 3rkd 

m V^kd/ V5kd/ 


(3.3) 


where Tq is the present time CMB temperature and qq and ^kd are the effective degrees of 
freedom of the cosmic plasma. Using Tq = 2.725 K = 2.35 x 10“^^ GeV, we have 


= 8.28 X 10“^® 



100 GeV 


m 


^ 20 MeV ^ 


(3.4) 


where we put in evidence the relevant numbers for the DM particle candidate, the neutralino. 
Note that the effective degrees of freedom qq and i^kd take into account the entropy conserva¬ 
tion when the cosmic plasma is heated by annihilation of particles since kinetic deconpling. 
Equation (3.4) shows that if a constraint is imposed on /3^ this is equivalent to say that we 
are imposing a constraint on the product mTkd, which represents a “finger-print” of the early 
interactions of DM particles with the cosmic plasma. 

We plot respectively in figs. 1 and 2 the CMB TT power spectrum and the DM power 
spectrum for /3^ = 10“^"^ and = 10“^® at 2 ; = 0. The standard CDM case (/3^ = 0) is 
also displayed for reference. In the lower panel of fig. 1 residuals with respect to our adopted 
reference (the standard ACDM model) are shown. With only a tiny correction (0.4 over 1000, 
for the = 10“^"^) the CMB spectrum does not effectively constrain the velocity dispersion 
of DM particles. In fig. 2 we compare the DM power spectra calculated for = 10~^® 
and 0^ = 10“^^ with those derived from WDM particles with masses m^dm = 1-4 keV and 
mw dm = 3.3 keV. For the plots of the WDM power spectrum we used the fitting formula 
provided by [20]. The vertical solid magenta and cyan lines (for /3^ = 10“^^ and (3“^ = 10“^®, 
respectively) represent the scales k = 21.7 and k = 68.5 h Mpc“^, respectively, for which 
0^ = at horizon crossing, thus constituting the limit of validity of our equations. 

The masses of WDM particles were chosen from Lyman-a forest constraints [20, 21]. 
Then we computed the values for which produce similar cut-offs. The plots in fig. 2 
indicate that velocity dispersion effects in GDM may mimic the cut-off due to free-streaming, 
typical of WDM particles. En passant, it is worth noting that, albeit such a cut-off is still 
consistent with constraints based on Lyman-a forest data, for masses larger than 3.3 keV, 
WDM appears no better than GDM in solving the small scale GDM anomalies [25]. 
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Figure 1. Upper panel. The CMB TT power spectrum computed for /3^ = 10 and /3^ = 10 
Lower panel. Residuals with respect to the reference ACDM model, represented with a dotted line. 



Figure 2. The DM power spectra at z = 0 computed for = 10“^® (cyan dashed line) and 
= 10“^"^ (magenta dashed line) compared with the corresponding ones for WDM particles of mass 
3.3 keV (cyan dotted line) and 1.4 keV (magenta dotted line). For the latter we used the fitting 
formula of [20]. The vertical magenta and cyan solid lines determine the limiting scales, defined as 
those for which at horizon crossing. 


It is important to emphasize that although in the considered case WDM and CDM 
produce similar cut-offs in the linear power spectrum, different mechanisms are involved: 
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in the case of WDM an exponential cut-off is produced in the distribution function due to 
free-streaming [41, 42] while velocity dispersion effects of CDM particles act as a pressure 
perturbation [18]. 

In order to mimic the cut-off in the linear power spectrum produced by 3.3 keV WDM 
particles it is required that = 10“^^. This condition cannot be satisfied by non-relativistic 
DM thermally produced. 

Therefore, let’s consider the possibility that DM particles were produced non-thermally. 
In this case, fixes the instant when these particles become non-relativistic. Indeed, if 
we assume that the transition to the non-relativistic regime occurs when cJq ~ 0 .1, then 
/3^ = 10“^® implies that the transition happens when Zc ~ 10^. 

Although this requires some fine-tuning in the production mechanism in order that pri¬ 
mordial nucleosynthesis be not perturbed, this possibility may offer an issue for the problem 
of the excess of satellites. 

If DM has a non-thermal origin able to explain, for instance, a parameter /3^ = 10“^®, 
we have seen that the non-relativistic phase should begin around Zc ~ 10^. This is similar 
to the Late Forming Dark Matter scenario (LFDM) proposed in ref. [43], see also [44]. In 
the LFDM picture a scalar field trapped in a metastable vacuum begins to oscillate and 
suffers a transition in which its equation of state passes from “radiation” {w ~ 1/3) to 
“cold matter” {w ~ 0). The matter power spectrum in this model typically shows a sharp 
break at small scales, below which power is suppressed [43, 45]. Cosmological simulations 
adopting this scenario were performed in [45], who have assumed that the transition occurred 
at zt = 1.5 X 10®, comparable to the redshift required to have /3^ = 10~^®. Comparing our 
fig. 2 with fig. 1 of reference [45], in which is shown the matter power spectrum resulting from 
the LFDM model as well as those resulting from WDM and ACDM models, it is interesting 
to notice that the present model, including velocity dispersion effects, produces a power 
spectrum similar to the LFDM picture. 

In fig. 3 is plotted the density of halos above a given mass scale M as a function of 
the mass for z = 0. The integral mass spectrum was computed for = 0 (ACDM model), 
/3^ = 10“^® and /3^ = 10“^^. For comparison, it is also plotted the integral mass spectrum 
for WDM models having particles with masses respectively equal to 1.4 keV and 3.3 keV. 
Note that for the case = 10“^® there is a slight deficiency of halos for masses less than 
2 X 10® Mq, comparable to the WDM model with 3.3 keV particles. In fig. 4 it is shown, for 
the same aforementioned parameters, the redshift at which the mass variance is equal to the 
unity as functions of the scale. A top-hat filter normalized to the Planck 2015 result, that is 
as = 0.829 ± 0.014 [46] was used in the computations. 

3.1 The mass of dark matter particles in the case /3 = 10“^® 

We have seen that = 10“^® can be explained only if dark matter particles have a non- 
thermal origin otherwise the predicted cosmic abundance would be inconsistent with the 
observed one. In this case, the fiducial value that hxes the initial conditions for the dispersion 
velocity corresponds to the transition between relativistic and non-relativistic regimes or 
corresponds to a “late formation of (non-relativistic) dark matter particles”. In both cases, 
the fiducial value corresponds to a critical redshift of about ~ 10^ or to a temperature of 
about 2.4 keV, thus after the primordial nucleosynthesis. 

What is the mass of particles related to those possible scenarios? Just after the critical 
redshift the distribution function obeys the Einstein-Vlasov equation whose solution for a 
flat FRW model is of the form / = f{ap/po) with ap/po = constant. The constant po has the 
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Figure 3. Number count as function of the mass scale computed for /3^ = 10“^® (cyan dashed line) 
and /3^ = (magenta dashed line) compared with the corresponding ones for WDM particles of 

mass 3.3 keV (cyan dotted line) and 1.4 keV (magenta dotted line). 



Figure 4. Redshift at which the mass-variance attains unity as function of the scale of the fluctuation 
computed for = 10“^® (cyan dashed line) and 0^ = (magenta dashed line) compared with 

the corresponding ones for WDM particles of mass 3.3 keV (cyan dotted line) and 1.4 keV (magenta 
dotted line). 


dimension of a momentum and was introduced in order that the argument of the distribution 








function be dimensionless. In this case, the density of dark matter particles is given by 


n = 


1 /‘°° 

2 ^/ dpp^f(ap/p„) = ^I,. 


where we introduced the integral fnnction 


PC 

Ik= dx x^f{x) . 

Jo 

The dimensionless velocity dispersion, in the non-relativistic case, is given by 


a^ = 


I roo 

2aX ■ipp"^/{ap/po) = 


m?a? I2 


Thus, the parameter fJP is simply 


o2 2 2 Po -^4 
p = (T^a^ = -^-r 
I 2 


(3.5) 


(3.6) 


(3.7) 


(3.8) 


In order to estimate po we require that these particles have the adequate cosmic abundance, 
in other words they have to satisfy 
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Eliminating poj one obtains: 


/3' 


= 3.45 X 10"^® 


h 

7 - 5/3 

-'2 





(3.9) 


(3.10) 


This resnlt is general, being the details of the particle distribution function encoded in the 

/ Q 

ratio 1 which involves integrals depending on the distribution function. 


4 Conclusions 

In this paper we report an study on velocity dispersion effects on the evolution of small 
perturbations of DM density. We considered the perturbed Vlasov-Einstein equation for 
DM an computed its momenta up to the second one, truncating the hierarchy based on the 
“quasi-cold” condition, for which we neglect terms 0{j?jrr?'). We implement this set of 
equations in GAME thereby computing the CMB TT power spectrum as well as the matter 
power spectrum. 

For DM candidates issued from extensions of the standard model like the neutralino, 
no signihcant velocity dispersion effects are expected either in the CMB TT power spectrum 
or in the matter power spectrum. Interactions with leptons keep these particles coupled 
thermally until temperatures of the order of 10-20 MeV. Consequently the expected value 
of the parameter is around 10“25^ several orders of magnitude smaller than the values 
required to affect the CMB TT power spectrum or the linear power spectrum. 

However, if DM particles are produce lately (zc ~ 10^) and non-thermally, higher values 
of the parameter are possible. In this case, velocity dispersion effects introduce a cut-off in 
the linear matter power spectrum at small scales, similar to that produced by WDM having 
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particles of mass around 3.3 keV. It is important to emphasize that the physical nature of 
the damping effects due to WDM and velocity dispersion of CDM is not the same: the 
former is related to “free-streaming” while the latter mimics pressure effects, see ref. [18] 
for a discussion on this point. The present model is comparable to the Late Forming dark 
matter scenario, in which non- relativistic matter appear as consequence of the decay of a 
scalar field. 
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